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1 Introduction and summary 



Understanding the spin structure of the nucleon is a central topic of present nuclear and particle physics 
activities, for a review see Of particular interest are certain sum rules which connect information at all 
energy scales, like e.g. the Gerasimov-Drell-Hearn (GDH) sum rule and its generalization to finite photon 
virtuality or the Burkhardt-Cottingham (BC) sum rule. Such sum rules are interesting from the theoretical 
point of view because they constitute moments of the sought after nucleon spin structure functions g\ and g-i- 
On the experimental side challenging new meson production experiments using real or virtual photons play an 
important role since only recently it has become possible to work with polarized beams and polarized targets, 
thus offering the possibility of mapping out the nucleons' spin structure encoded in these two functions, which 
can be formulated on a purely partonic (high energy regime) or hadronic level (low energy regime). In both these 
extreme cases, systematic and controlled theoretical calculations can be performed. The region of intermediate 
momentum transfer is accessible using quark/resonance models or can be investigated using dispersion relations. 
In fact, one of the final goals of many of these investigations is to obtain an understanding of how in QCD this 
transition from the non-perturbative to the perturbative regime takes place, guided by the precise experimental 
mapping of spin-dependent observables from low momentum transfer to the multi-GeV region, as it is one of 
the main thrusts of the research carried out e.g. at Jefferson Laboratory. 

Here we focus on a theoretical investigation of the nucleon's spin structure in the non-perturbative regime 
of QCD. We are utilizing chiral perturbation theory (CHPT) to analyze the structure of the nucleon at low 
energies, based on the spontaneous and explicit chiral symmetry breaking QCD is supposed to undergo (for 
a general review, see e.g. Q). By now it is well established that the pion cloud plays an important role in 
understanding the nucleons properties in the non-perturbative regime of QCD, and many processes have been 
analyzed using chiral perturbation theory. Some recent work in various versions/extensions of baryon CHPT 
pertinent to the topics discussed here can be found e.g. in Refs. || ||, ||, |[ [?], ^] . In the letter || we had presented 
a novel analysis of chiral loop effects in the generalized Gerasimov-Drell-Hearn sum rule in the framework of 
a Lorentz-invariant formulation of baryon chiral perturbation theory |l0| . We performed a complete one-loop 
calculation for the spin-dependent proton and neutron structure functions S\p,(y = 0,Q 2 ) of forward double 
virtual Compton scattering (V 2 CS), where v is the energy transfer and Q 2 the negative of the photon virtuality 
(momentum transfer squared). Combining the methods of CHPT and analyticity, one can deduce all desired 
spin-related sum rules, respectively all moments of the nucleon's spin structure functions, from these two basic 
building blocks for Q 2 < 0.3 GeV 2 . Results where thus presented for the V 2 CS structure functions S±' n (0, Q 2 ) 
(where the bar indicates the subtraction of the elastic intermediate state) and for the first moment of the sum 
rule for the spin-dependent structure function g\ (for precise definitions, see below). Also, a comparison to 
the previously obtained so called "heavy baryon" results, which are implicitly contained in our calculation and 
can be extracted by performing a 1/ (nucleon mass) expansion of our relativistic results, was given. Here, in 
the light of the sometimes rather slowly converging 1/m expansion of the heavy baryon framework in the spin 
sector ||, we do not only give a detailed exposition of how the relativistic results have been obtained, but 
also give novel results for the structure functions ^'"(O, Q 2 ), for the generalization of the DHG and BC sum 
rules as well as for the forward and longitudinal-transverse spin polarizabilities. These latter two quantities 
parameterize the forward spin-dependent (virtual) Compton scattering amplitude at low energies and encode 
important information about higher moments of the nucleon's spin structure functions. 

The pertinent results of this investigation can be summarized as follows: 

1) We have studied spin-dependent doubly virtual Compton scattering off nucleons to fourth order (one 
loop) in a Lorentz-invariant formulation of baryon chiral perturbation theory. At this order in the chiral 
expansion, no unknown low-energy constants appear and thus we obtain parameter-free predictions for 
the spin structure functions Si^iy, Q ), where the bar signals the subtraction of the elastic intermediate 
state. These two V 2 CS structure functions serve as the fundamental building blocks from which all other 
results can be derived. 

2) We have studied the corresponding sum rules (the generalization of the Gerasimov-Drell Hcarn sum rule 
to finite photon virtuality and the Burkhardt-Cottingham sum rule) and their moments, which can be 
obtained by expanding the structure functions at low energies is, cf. Eqs. ( ^26| . ^27lj2~30| ). 
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3) We have estimated the contribution of A(1232) Born graphs and from vector meson contributions as 



described in Section 3.3. The delta plays an important role in the nucleon spin sector and thus an 



extension of our work including nA loops in a systematic manner is called for. 

4) Results for S?' n (0, Q 2 ) and the first moment Y\' n {Q 2 ) as well as for T\~ n {Q 2 ) were already given and 
discussed in M . We note that the prediction for T\~ n (Q 2 ) is free of spin-3/2 effects and in rough agreement 
with preliminary data from Jefferson Laboratory plj . Predictions for the structure functions Sf' n (0, Q 2 ) 
and for the moment r^(Q 2 ) including resonance contributions, as well as for the chiral expansion of the 
second moments s[ 2 ^ p ' n (0, Q 2 ) are shown in Fig. [|, Fig. [| and Fig. [|, respectively. 

5) The one-loop results for the second structure function are shown in Fig. ^ for the proton and the neutron. 
They come out sizeably different from the heavy baryon result for the same reasons as already discussed 
in ||] in connection with the first structure function. The delta Born and vector meson contributions have 
also been evaluated, see Fig. 0. 

6) The chiral prediction for a linear combination of first moments of g\ and g2 commonly denoted as I^(Q 2 ) 
supplemented with resonance contributions (mostly from the delta) agrees with recent data from Jefferson 
Laboratory |l2| for virtualities Q 2 < 0.3 GeV 2 , see Fig. ||. However, if one includes a ANj vertex function 
as suggested by older electroproduction data at higher energies, one is not able to describe the curvature 
of the data as the virtuality increases, see also Fig. ||. 

7) The Burkhardt-Cottingham sum rule is fulfilled within the Lorentz-invariant approach used here. How- 
ever, truncating the product of form factors at fourth order does not give a very accurate representation of 
the full result, which can be obtained using the form factors from dispersion theory or infrared regularized 
baryon CHPT supplemented by vector mesons |T^| . 

8) We have calculated the forward spin polarizability 70 (Q 2 ) (which corresponds to a linear combination 
of the second moments of g\ and g-i) at zero and finite photon virtuality to 0(p A ) in relativistic baryon 
CHPT for the proton and the neutron. To obtain a reliable prediction for the chiral corrections to these 
spin polarizabilities at Q 2 = in the heavy baryon approach, we find that one has to account (at least) for 
terms of 0(M^) : i.e. go to N 4 LO in the calculation. However, we note that even our relativistic one-loop 
result, when supplanted with resonance model delta Born term and vector meson contributions, still does 
not agree with the data |14fl . We have also isolated the reason for the bad convergence, it rests entirely 
on the slow convergence of the pole (Born) graphs. 

9) We have calculated the longitudinal-transverse spin polarizability Sq(Q 2 ) (which corresponds to a linear 
combination of the second moment of g\ and the first moment of g^) at zero and finite photon virtuality 
for proton and neutron. When performing a chiral expansion on the full relativistic results, including 
terms up-to-and-including G{M 7! ) 1 one finds a good convergence for these quantities at the photon point 
in the 1/m expansion. This observation and the absence of large delta (Born-term) contributions leads to 
the hope that 6o(Q 2 ) can serve as a testing ground of the chiral dynamics of QCD which might be better 
suited than the forward spin polarizabilities of the nucleon. 

f 0) We have calculated the fifth order contributions from the pole (Born) graphs for a particular combination 
of the structure functions of the neutron shown in Fig. [H^]. The complete fifth order result for the neutron 
BC sum rule is in good agreement with the product of the form factors for Q 2 < 0.f5GeV 2 , cf. Fig. [l4| 
The partial fifth order contribution to the neutron's spin polarizability is non-negligible and negative, 
indicating that these terms might be the origin for explaining the discrepancy to the data in case of the 
proton. 



The material in this paper is organized as follows. Section |2| contains the basic formalism for spin-dependent 
doubly virtual Compton scattering off nucleons, its relation to pion electroproduction and the derivation of the 
pertinent sum rules and their moments. In section |[ the theory underlying our calculations is exposed together 
with the chiral expansion of the spin-dependent structure functions. We also show how we estimate various 
resonance contributions. The results are presented and discussed in section [|. Some technicalities are relegated 
to the appendix. 
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2 Basic formalism 



This section contains the basic formalism on spin-dependent doubly virtual Compton scattering off nucleons, its 
relation to inelastic electroproduction and the derivation of the corresponding sum rules. The reader familiar 
with these concepts might directly proceed to section ||. 

2.1 Doubly virtual Compton scattering 

We consider spin-dependent doubly virtual Compton scattering (V 2 CS) off nucleons (neutrons or protons) in 
forward direction, that is the reaction 

7 *(«, e) + N(p, s) - j*(q, e') + N(p, a') , (2.1) 

with q (jp) the virtual photon (nucleon) four-momentum, s (s') the nucleon spin (polarization) and e (e') the 
polarization four-vector of the incoming (outgoing) photon. It is common to express the spin amplitude of 
V 2 CS in terms of two structure functions, called S\(v,Q 2 ) and S%(v, Q 2 ), via 

TlH (p, q , s) = -ij± e^q a | S/3 S^Q 2 ) + [p ■ q Sp - S ■ q p ] } , (2.2) 

where s M denotes a spin-polarization four-vector, m is the nucleon mass, e^ af3 the totally antisymmetric Levi- 
Civita tensor with e 0123 = 1, v = p ■ q/m the energy transfer and Q 2 = —q 2 > the (negative of the) photon 
virtuality. N is a convention dependent normalization factor. Following Ref.|j| we set 5*1(0,0) = —e 2 n 2 /m 2 . 
With the spinor normalization uu = 1 typically used in baryon CHPT, we find N = 2m, which we will utilize 
throughout this work. Note that while S\(y, Q 2 ) is even under crossing v <-> — v, the structure function 5*2 (y, Q 2 ) 
is odd. 

In the rest frame of the nucleon the forward virtual Compton tensor in the Coulomb gauge eo = e' = simplifies 
to 

rife Q 2 ) = \ e ijkX ^ {1/^1(1/, Q 2 ) + [v 2 o k -a-qq k ] X , (2.3) 

where v now corresponds to the photon energy in the laboratory system, the (k — 1,2,3) are the Pauli 
spin-matrices and x is a two-component spinor. Utilizing the identity 

a ■ qe' ■ (e x q) = (u 2 + Q 2 ) [a ■ (e" x e) - {a- (e' x q) e ■ q - a ■ (e x q) ? ■ q)] (2.4) 

we finally arrive at the V 2 CS forward matrix element 

e'-T.e| rest = ^x 1 {™-(?' x?)[mvS 1 (v,Q 2 )-Q 2 S 2 (v,Q 2 )] 

+i [a ■ (e' xq)e-q-a- (ex q) t ■ q] {v 2 + Q 2 ) S 2 (v, Q 2 )} \ , (2.5) 

in the form which is best suited for our calculations in the framework of relativistic baryon chiral perturbation 
theory. The Compton amplitudes S\.2{v, Q 2 ) are amenable to a chiral expansion as will be detailed in section|3[ 
We remark that in what follows, we will mostly be concerned with the reduced amplitudes 

Si(v, Q 2 ) = Si(u, Q 2 ) - Sf(v, Q 2 ) , (2.6) 

i.e. the Compton amplitudes with the contribution from the elastic intermediate state subtracted. More 
precisely, these are the contributions from the single nucleon exchange (pole) terms with the corresponding 
vertices given in terms of the electromagnetic form factors. Only the non-pole parts of the corresponding 
diagrams contribute to the nucleon spin structure as discussed in more detail below. We note that for the 
calculations to follow we mostly work at v = (and therefore often do not make explicit the energy dependence 
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of certain quantities). In fact, for the calculations we make use of crossing symmetry, which allows to write the 



relativistic V 2 CS forward matrix element in terms of two invariant functions A(v, Q 2 ) and B(is, Q 2 ) &t# 5 : 
e'-T-e\ rcst = X t{#.(^ xg ) [A(s,Q 2 )-A(2m 2 -2Q 2 -s,Q 2 )] 

+ i[ff. (e' xq)e-q-a- (ex q) e" ■ q] \q \ 2 [B(s, Q 2 ) - B(2m 2 - 2Q 2 - s, Q 2 )} } X , 

(2.7) 

with s = (p + q) 2 the standard Mandelstam variable, which is related to v via v = (s — m 2 + Q 2 )/2m. Written 
in term of v, crossing symmetry means that in this equation the difference of a quantity taken at v with the 
same quantity taken at —v appears. Since we consider quantities at v = i.e. at s = m 2 + q 2 , all of these can 
be obtained as derivatives of order n > 1 of the functions A and B. More precisely, we give two examples to 
illustrate this method. 



afroygiim^ '^ = 4 dA (^ q 2 ) 

m v— >o v dv 
Q 2 hm M0,Q 2 ) A dB{y, Q 2 ) 



(2.8) 

!/=0 

(2.9) 



li/=0 ' 



iy^o 1/ dv 

where j/5 2 (^, Q 2 ) = vS 2 {v, Q 2 ) - vS 2 {v, Q 2 ) 
2.2 Relation to inelastic electroproduction off nucleons 

Because of unitarity, there is a basic connection between spin structure functions in V 2 CS and the ones probed 
in inelastic electroproduction experiments. This is simply related to the fact that the imaginary part of the 
Compton tensor is given in terms of nucleon plus meson states, the lowest one being the pion-nucleon state. 
The differential cross section for exclusive electroproduction can be expressed in terms of four virtual pho- 
toabsorption cross sections ot{v,Q 2 ), gl{v,Q 2 ), c'ltI^'Q 2 ) ano - ®tt( v iQ 2 )- The transverse, <Tt, and the 
transverse-transverse, cr' TTl cross sections can be related to the total photoabsorption cross sections <J\i<i and 
CT3/2 corresponding to the excitation of intermediate states with spin projection 3/2 and 1/2, respectively: 

°t = 2 ( cr 3/2 + 01/2) ) 

a TT = 2 ( CT 3/2 - CT1/2) ■ ( 2 -10) 

The four virtual photoabsorption cross sections are related to the quark structure functions Ft, F2, Gi and G2 
(that depend on v and Q 2 ) measured in (un)polarized deep inelastic electron scattering off nucleons via: 

Air 2 rv 

*t(v,Q 2 ) = —^F^Q 2 ), (2.11) 
mK 

_ fAffifilOTp^.flfil^)!, (2 . 12) 



K 

47r z a ,„ „ v 



v'lA^Q 2 ) = r^7^(Gi(^,Q 2 ) + -G 2 (^,Q 2 )) , (2.13) 

K m 

v'ttM 2 ) = -^v(Gi(v,Q 2 )~-j 2 G 2 (v,Q 2 )) , (2.14) 
K m 

where v = E — E' is the virtual photon energy in the lab frame, a = e 2 /An = 1/137.036 the fine structure 
constant, 7 = Q/v and K a flux factor. For more details on the electroproduction formalism, see e.g. Ref . Q . 

2.3 Sum rules and moments 

Let us first consider the real photon case with Q 2 = 0. At low energies one expands the coefficient 

2mj(i/)= Ynn mvS l {v 1 Q 2 )-Q 2 S 2 {v,Q 2 ) (2.15) 



# 5 A similar method was used e.g. in Rcfs. Jl5|, ^ \^ . 
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in front of the spin-dependent operator a • (e ' x e*) in powers of the energy transfer, 

2 2 

giy) = -y^v + 10 + 0{v 5 ) , (2.16) 

where the first term proportional to the nucleon anomalous magnetic moment squared (~ k 2 ) is given by 
the venerable low-energy theorem (LET) of Low, Gell-Mann and Goldberger |L7| and 70 is the forward spin 
polarizability. The imaginary part of g{v) can be related via the optical theorem to 0\j 2 and 03/21 

Im g{v) = V - {<j 1/2 {v, Q 2 = 0) - a 3/2 (v, Q 2 = 0)) (2.17) 



Writing an unsubtracted dispersion relation for g(y) one obtains using Eq.(2.16) and (2.17) the celebrated 
Drell-Hearn-Gerasimov sum rule Jl8fl (we give here the generic form without specifying whether we mean the 
neutron or the proton): 

2 r K 2 ol k 2 f dv r n . „ . 

— h/2-CT 1/2 ] . (2.18) 

The value of the LHS of the sum rule is —204 [ih for the proton and —232 [ib for the neutron. This sum rule is 
presently under active experimental investigation at MAMI, ELSA, GRAAL, CEBAF and Spring-8. For first 
results, see Jl9[ . 

The DHG sum rule has been generalized to finite Q 2 [^o|. It is most natural to connect such a generalization 
of the sum rule to the Bjorken sum rule in deep inelastic scattering. One thus defines 



2m 2 



HQ )=qT J 9i(x,Cf)dx, (2.19) 

with x corresponding to the pion production threshold. Here, G\ — gi/{mv) can be related to CT3/2 and ai/ 2 
via Eqs.( 2.l6] ,2.14) and x = Q 2 /2mv is the standard scaling variable. At zero virtuality, one recovers the DHG 
sum rule, i.e. ii(0) = — k 2 /4. Similarly, there exists a sum rule for the second spin structure function, g 2 , the 
so-called Burkhardt-Cottingham (BC) sum rule 

x 



2m 2 

Jo 

1G M (Q 2 )-G B (Q 2 )^ , n , 1 m2 . 



HQ 2 ) = I 9 2(x,Q 2 )dx (2.20) 



G M {Q 2 ) = -F 2 {Q 2 ) (F^Q 2 ) + F 2 (Q 2 )) , (2.21) 

\^ + 4m 2 I 

where Ge(Q 2 ) and Gm(Q 2 ) are the electric and the magnetic nucleon Sachs form factor, respectively, Fi(Q 2 ) 
and F 2 (Q 2 ) are the Dirac and the Pauli form factor, in order, and G 2 — gijv 2 . We note that in a Lorentz- 
invariant calculation as performed here, the Dirac and Pauli form factors appear naturally, while in the HBCHPT 
approach the Sachs form factors are the more appropriate quantities f22|| . This sum rule assumes that a' LT 
vanishes faster than 1/u at large v. Other sum rules have been defined in the literature, see e.g. || [l6|| : 



2m 2 



x 



Ia(Q 2 ) = -@rj (9i(x,Q 2 )-Yg2(x,Q"))dx , (2.22) 
2m 2 r x ° 1 

Ic(Q 2 ) - j—- } (g 1 (x,Q 2 )-j 2 g 2 (x. / Q 2 ))dx. (2.23) 

All these sum rules can be written in terms of the structure functions S\(i/, Q 2 ), S 2 (v, Q 2 ) of V 2 CS using the 
following dispersion relations : 



S^Q 2 ) = 4e 2 / 
Jc 



dv'v'Gxiy' \Q 2 ) 



Q2 j 2m V — V 

JQ 2 /2m V* - V 1 



G 



where use has been made of the optical theorem, 

Im S t (v, Q 2 ) = 2tt Gi(v, Q 2 ) , (i = 1,2) . (2.25) 

Expanding the structure functions at low energies v, that is around v = (note that we now consider the 
structure functions with the elastic intermediate state subtracted), 

§i(v, Q 2 ) = S[ 0) (0,Q 2 ) + u 2 S[ 2) (0,Q 2 ) + u 4 S[ 4) (0,Q 2 ) + ... , 
vS 2 {v,Q 2 ) = 5^ 1) (0,Q 2 ) + ^ 2 ^ 1) (0,Q 2 ) + ^ 4 ^ 3) (0,Q 2 )+^^ 5) (0,Q 2 ) + ... , (2.26) 
where the ellipsis denote terms of higher order in i/, one obtains the following set of sum rules for S%(0, Q 2 )- 

s<W) . ..r^i.^^, (2 . 27) 

Jvo v m 

S?\0,Q 2 ) = r ^XiP (< = 1.2.3,.-.). (2-28) 

noting that 5i°^(0, Q 2 ) is directly related to Ii(Q 2 ). Note that the often used first moment Ti(Q 2 ) is related to 
h(Q 2 ) via 

T 1 (Q 2 ) = ^h(Q 2 ) . (2.29) 



Likewise, 5^ (0, Q 2 ) satisfies: 
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SfV&Q 2 ) = 4e 2 / di/G 2 {i/,Q 2 ) = ~^h{Q 2 ) , (2.30) 
#(0,Q 2 ) = 4e 2 r « = 1,3,5,...), (2.31) 

where 1 ^(0,Q 2 ) is related to I 2 - In order to calculate 7c it is necessary to define a new quantity Sr similar 
to Si — Q 2 S2/v, but having only right hand cuts: 



S R { V , Q z ) = 4e^ / &/ ±v y v y/ , (2.32) 
so that the following relations hold: 

Ia(Q 2 ) = ^(^ 0) (0 > O a )-O a ^ 1) (Q,Q 2 )) , (2.33) 

9 

777 _ 

J C (Q 2 ) = — S R (Q 2 /2m,Q 2 ). (2.34) 



In Eq.(2.1£) the forward spin polarizability appears j22| ], which is defined in terms of the transverse-transverse 
photoabsorption cross section as 

70 = i/ ~ ^V^M) - ^3/2(^,0)) . (2.35) 

Note that x is trivially zero at vanishing virtuality in this equation so that one recovers the usual dispersion 
relation for 70. One can define a similar quantity involving a' LT instead of <r' TT , the so-called longitudinal- 
transverse polarizability, denoted So in what follows, 

^ = i/$(l-x) Q lim(^. 1/2 (,,0 ) ). (2.36) 

These two quantities can be generalized to finite Q 2 . The quantity x was in fact introduced in order to get 
straightforward relations between 7o(Q 2 ), <5o(Q 2 ) and the Si(iy,Q 2 ). One has: 

7o(Q 2 ) = ±(s{ 2 \o,Q 2 )-^si 3 \o,Q 2 )) , (2.37) 

57T V TO j 

1 

8^ 



MQ 2 ) = ^($ 2) (0,Q a ) + s£ 1) (0,Q 2 )) . (2.38) 
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3 Chiral expansion of the structure functions 



Our calculations are based on an effective chiral pion-nucleon Lagrangian in the presence of external sources 
(like e.g. photons) supplemented by a power counting in terms of quark (meson) masses and small external 
momenta. Its generic form consists of a string of terms with increasing chiral dimension, 

C « - + £ (2) + £ (3) 4- £ (4) + £ (2) + £ (4) + (31) 

The superscript denotes the power in the genuine small parameter q (denoting pion masses and/or external 
momenta). A complete one-loop (fourth order) calculation must include all tree level graphs with insertions 
from all terms given in E q. (|3.l[ ) and loop graphs with at most one insertion from C^ N . The complete Lagrangian 
to this order is given in |23fl. We note that for the case under consideration the only appearing dimension two 
low-energy constants (LECs), called cq and cj 111, can be fixed from the anomalous magnetic moment of 
the proton and of the neutron, respectively. As discussed below, there are no contributions from for 
the observables considered here. Besides these chiral corrections to V 2 CS, we also estimate various resonance 
contributions to the observables, however, we do not have an equally systematic framework at our disposal for 
that task. So some future work should be undertaken to sharpen this part of our calculation by e.g. extending 
the small scale expansion of p5| to Lorentz-invariant formulation. 



3.1 Essentials of Lorentz-invariant baryon CHPT 

Baryon chiral perturbation theory is complicated by the fact that the nuclcon mass does not vanish in the chiral 
limit and thus introduces a new mass scale apart from the one set by the quark masses. Therefore, any power 
of the quark masses can be generated by chiral loops in the nucleon (baryon) case, spoiling the one-to-one 
correspondence between the loop expansion and the one in the small parameter q. One method to overcome 
this is the heavy mass expansion (called heavy baryon chiral perturbation theory, for short HBCHPT) where 
the nucleon mass is transformed from the propagator into a string of vertices with increasing powers of 1/m. 
Then, a consistent power counting emerges. However, this method has the disadvantage that certain types 
of diagrams are at odds with strictures from analyticity. The best example is the so-called triangle graph, 
which enters e.g. the scalar form factor or the isovector electromagnetic form factors of the nucleon. In a 
fully relativistic treatment, such constraints from analyticity are automatically fulfilled. It was recently argued 
in [£4[ that relativistic one-loop integrals can be separated into "soft" and "hard" parts. While for the former 
the power counting as in HBCHPT applies, the contributions from the latter can be absorbed in certain LECs. 
In this way, one can combine the advantages of both methods. A more formal and rigorous implementation of 
such a program was given in fl(i|| . The underlying method is called "infrared regularization" . Any dimensionally 
regularized one-loop integral H is split into an infrared singular (called I) and a regular part (called R) by a 
particular choice of Feynman parameterization, 

H = I + R. (3.2) 

Consider first the regular part. If one chirally expands the contributions to R, one generates polynomials in 
momenta and quark masses. Consequently, to any order, R can be absorbed in the LECs of the effective 
Lagrangian. On the other hand, the infrared (IR) singular part has the same analytical properties as the full 
integral H in the low-energy region and its chiral expansion leads to the non-trivial momentum and quark-mass 
dependences of CHPT, like e.g. the chiral logs or fractional powers of the quark masses. For a typical one- 
loop integral (like e.g. the nucleon self-energy E) this splitting can be achieved in the following way (omitting 
prefactors) 

, d d k 1 f 1 , f d d k 1 

^ = / ~, r-?—. — = I dz 



(27r) d AB J Q J (2ir) d [(1- z)A + zB} 2 

oo /-cxn r d d k 



with A = M 2 — k 2 — ie, B — m 2 — (p — k) 2 — ie, e — > + , M w the pion mass and d the number of space- 
time dimensions. Any general one-loop diagram with arbitrary many insertions from external sources can be 



dZ l (2nyi(l-z)A + zB) 2 =I + R > (3 - 3) 



S 



brought into this form by combining the propagators to a single pion and a single nucleon propagator. It was 
also shown that this procedure leads to a unique, i.e. process-independent result, in accordance with the chiral 
Ward identities of QCD jl0| . Consequently, the transition from any one-loop graph H to its IR singular piece 
/ defines a symmetry-preserving regularization. For more details, the reader is referred to [fiof . 

There is one special feature to IR regularized baryon CHPT that we have to discuss. Here, we are interested 
in extending the range of photon virtualities Q 2 to values that eventually allow a matching to the perturbative 
QCD calculations. However, there is a generic limit to such an extension. A priori, infrared regularization defines 
a low energy effective field theory, but due to the resummation of kinetic energy insertions and similar effects, 
one might hope to extend this framework to higher virtualities than the HBCHPT scheme (this expectation 
was borne out explicitely in the calculation of the neutron charge form factor, see [p^|). Consider now a specific 
one-loop function in d dimensions 



i J (2tt 



1 f d d l I 



i) d {Ml - I 2 ) (m 2 -(p + q- I) 2 ) (m 2 - (p - I) 2 ) 
with (p + q) 2 = s and p 2 = m 2 . In infrared regularization, this integral is written as 



(3.4) 



r ( S ,Q 2 )=r(3-^ / dy I dzzC%- 3 (y,z) , (3.5) 



2 



o Jo 



with C a polynom in y and Z, 



and 



C = Az 2 + Bz + /i 2 , 
A = l + 2yfl[i + y f i 2 -y{l-y)T , 

B = -2yfi/i-(l + 2;)/i 2 , (3.6) 

a - m 2 - Ml M„ q 2 Q 2 
Q = o — /Tr ' ^ = ' T = ~~ 2 = 2 ■ ( 3 - 7 ) 

The T-function in Eq^^) contains the pole as d ~ + 4 which has to be dealt with in standard fashion. We also 
work with the standard value of the scale of dimensional regularization setting X — m. All quantities considered 
here contain derivatives with respect to s of such loop functions at s — m 2 + q 2 < m 2 (see the discussion at the 



end of Section 2.1). At that point, 

A = l + yr-y(l-y)T= 1 + y 2 t , (3.8) 

so that A = for t = —1/y 2 . Consequently, r (s, Q 2 ) does not converge for |t| > 1. Thus, the photon virtuality 
is bounded by the nucleon mass. In fact, for higher derivatives the influence of this singularity is felt at smaller 



values of Q 2 so that we restrict the range of photon virtualities to Q, 2 nax = 0.3 .. . 0.5 GeV 2 depending on the 
observable one looks at. 



3.2 Structure functions 

Within this approach, we have calculated the reduced structure functions s[ P 2 n \o,Q 2 ), generically called S in 
this section. The chiral expansion of S takes the form 

S = S tTCC + S loop . (3.9) 

In our case, the tree level contribution stems from the remainder of the Born graphs which lead to the following 
amplitudes (we omit writing down the crossed term contributions) 

T Bom = CWJ_ + n j (31Q) 

s — m 2 

with s = (p + q) 2 , C(Q 2 ) can be expressed in terms of the nucleon electromagnetic form factors and TZ denotes 
the non-pole (polynomial) remainder from the Born diagrams. Only this latter contribution survives the 
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subtraction of the contribution from the elastic intermediate state (nucleon pole term) . Note also that there is 
no local contribution from C^ N . This can be understood easily. For constructing a gauge-invariant operator 
corresponding to a two-photon observable, one needs two factors of the electromagnetic field strength tensor 
.F^j,, which leads to terms of fourth order. However, such terms are always proportional to spin-independent 
terms like e.g. e' • e, i.e. they correspond to spin-independent Compton scattering (see e.g. [^6]). To generate 
an operator with the additional factor containing the nucleon spin matrices, one needs at least one extra power 
in small momenta, i.e. such terms can only start at fifth order in the effective Lagrangian. 

The pertinent one -loop diagrams are shown in Figs. ^ and ^ for the third and fourth order contributions, 
respectively. Note that different to the heavy baryon approach, we do not need to consider dimension two 
insertions from the kinetic energy on the nucleon lines since this is done automatically in a Lorentz-invariant 
formulation. In fact, all such diagrams are counting as third order in the approach used here, the only genuine 
fourth order graphs are the ones with one insertion ~ cq^ (anomalous magnetic moment). Note that these 
are indeed the only low-energy constants which enter in a complete one loop calculation (whenever possible, 
we substitute these LECs by the appropriate combinations of the proton and neutron magnetic moments). 
Therefore, the Si 2(f) Q ) are already non-vanishing at third order, quite in contrast to the heavy baryon 
calculation. As a check, we remark that in the limit of a very large nucleon mass, one recovers the earlier 



However, as already found in the relativistic 



HBCHPT results of [§, |§ and also the recent results of 
calculation of the slope of the generalized GDH sum rule in [|3| , one can not give the results in closed analytical 
form. The explicit expressions for the one-loop contributions to the structure functions in terms of the pertinent 
loop functions are given in the appendix. Whenever possible, we will compare our results to the ones obtained 
within HBCHPT. It should already be noted here that we find a very strong dependence on the value of the 
nucleon mass, not untypical for the nucleon spin sector. Of course, there are also non-negligible resonance 
contributions to the Si(i>, Q 2 ). More precisely, the effects from the A are expected to largely cancel in Si(v, Q 2 ) 
H but are supposed to be more important in ^(f, Q 2 )- There are additional (smaller) vector meson and higher 
baryon resonance contributions. In the following section, we briefly discuss how to estimate the contributions 
from these states. 



3.3 Modeling resonance contributions 

It is well-known that the excitation of the A(1232) plays a significant role in the spin sector of the nucleon. 
One therefore would like to include the delta as a dynamical degree of freedom in the effective Lagrangian. 
At present, that can only be done systematically in the heavy baryon scheme treating the nucleon-delta mass 
splitting as an additional small parameter [p5[ . An effective field theory formulation for the relativistic pion- 
nucleon-delta system does not yet exist. Therefore, to get an estimate of the contribution of the A-resonance 
to the various spin structure functions discussed here, we calculate relativistic Born graphs. These are obtained 
using the well-known relativistic spin-3/2 propagator (with momentum P, running from v to \x): 

f + niA ( 1 „ „ 2P M P, Ptf - P v >f 



and the A — * A7 transition operator (for an in-coming photon with momentum k and an out-going nucleon 
with momentum p) : 

^75 (W- fcV+ (Y + -)(#- #07" V 3 (3.12) 

' 3 . (3.13) 



2m V 2 

ieg 2 



Am 2 



V ■ /c 7 5 ( e" - (X + \)<?f ) - e • p( - (X + l)#f 



2"' / "V 2 



Here, T is the ~ — ► | isospin transition operator and is the delta mass. The A — ► A7 transition depends on 
two off-shell parameters X, Y and two transition strengths g\ and g 2 , quantities which are not so well known. 
We stress that in an effective field theory approach such a dependence on off-shell parameters would be lumped 
in certain low-energy constants, however, in the absence of such a systematic analysis of delta effects we must 
resort to the procedure adopted here. Bounds on A, Y, g\ and g 2 have been given in Ref. |26|]: —0.8 < A, Y < 
0.4, 4 < g% < 5 and 4.5 < gi < 9.5. Here we constrain g\ to its large N c relation, g\ = 3(1 + n p — n n ) /2y/2 = 5.0 
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and use two sets of parameters, X = Y = 0.4, 172 = 4.5, and X = Y = —0.8, 52 = 9.5, respectively. We 
note that these bounds are very conservative, a more precise determination based on a combined reanalysis of 
spin-independent Compton scattering and pion electroproduction based on IR baryon CHPT would certainly 
lead to more stringent bounds. Furthermore, in order to take into account the fact that the A — ► Nj transition 
occurs at finite Q 2 we entertain the possibility of introducing a transition form factor GAN-yiQ 2 ), following p| 



cxp(-0.2Q 2 /GeV 2 ) 
(l + 1.41Q 2 /GeV 2 



as extracted from pion electroproduction in the A(1232)-resonance region, see Refs. p7|. Note, however, that 
the energies and virtualities of these data were sizeably larger than the range of values considered here. In a way, 
such a form factor subsumes some higher order effects, it any way only plays a significant role for virtualities 
Q 2 > 0.15 GeV 2 . We will come back to this issue when we discuss the result in the next section. Of course, 
there are also smaller contributions from higher baryon resonances, but we do not include them in this work. 

A less pronounced though important resonance contribution is related to the vector mesons. Again, a systematic 
EFT prescription how to include these degrees of freedom does not exist. We adopt here the procedure advocated 
in Ref. [^3). In the pion-nucleon EFT, any vector meson contribution is hidden in the values of the various LECs. 
However, the momentum dependence of the vector meson propagator is only build up slowly by adding terms 
of ever increasing chiral dimension. This can be done much more efficiently by including vector mesons in a 
chirally symmetric manner and retaining the corresponding dimension two counterterms, so that the LECs ce, C7 
are effectively replaced by 

F P M p 

C6 — > C6 + 9pNN Hp M2 _ - , 

9 p nn n p FpMp g^NN FuMu g<pNN n<p F^M^ , , 

C7 ^^Mp] + ^J^ + ^Mp- (3 - 15) 

Here, t is the invariant four-momentum squared and the remainders c@,c 7 account for physics not related 
to vector mesons. They have been determined from fitting the nucleons electromagnetic radii Jtt3|. All other 
parameters appearing in Eqs. ( 3.15| ) can be taken form the dispersion-theoretical analysis of Refs. ]28| . The values 



of these parameters as used here are M p = 770 MeV, F p = 152.5 MeV, g pNN = 4.0, n p = 6.1, M u = 780 MeV, 
F u = 45.7MeV, g uNN = 41.8, k u = -0.16, M = 1020 MeV, F = 79.0 MeV, g^, NN = -18.3 and k = -0.22. 



4 Results and discussion 

Before presenting results, we must fix parameters. We use gA — 1.267, F v = 93 MeV, M„ = 139.57 MeV, 
m — 938.9 MeV, k v = 3.706 and k s — —0.120. The latter two quantities include the LECs cq and C7. The 
pertinent delta and vector meson parameters were already given in the preceding section. We stress again that 
no unknown parameters enter our calculation so all the following results are true predictions to fourth order 
(more precisely, this refers to the chiral expansion of the various quantities). 



4.1 Complete fourth order results 



We first discuss the two structure functions, the various sum rules and the moments thereof. Results on the 
chiral expansion for 5f' ra (0, Q 2 ) and T\ ,n {Q 2 ) as well as for T\~ n (Q 2 ) were already given in ||, together with 
a comparison to the HBCHPT results. Here, we show in the Fig. || the one-loop predicti ons for 5f' n (0, Q 2 ) 
supplemented by the vector meson and delta contributions calculated as described in Section 3.3. The band as- 
sociated to the resonance contributions is due to the parameter variations of certain delta couplings as discussed 
before. Of course, to draw a final conclusion on the size of the resonance part one has to include it A loops in 
a systematic manner (for some attempts see ||, [| ||). No published data exist so far for a direct comparison 
with these predictions, for orientation we show the prediction for T\(Q 2 ) including resonance effects in Fig. 

# 6 We remark that in some terms of higher order we did not replace the combinations of the LECs eg and C7 by re p and K n 
which leads to slightly different results for the structure functions. 
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In Figs. |_and | we show the one-loop results for the first moment S^\Q 2 ) as defined in Section |2~3| and for 
(Q 2 /mis) ^2(1^, Q 2 )\v=o (this is the quantity which indeed enters the various sum rules) for the proton and the 
neutron. The results are similar to the ones presented for 5f : ™(0, Q 2 ) in ||. The fourth order corrections are not 
overly large for virtualities below Q 2 ~ 0.15 GeV 2 and are very different from the ones obtained in HBCHPT. 
These differences can be understood from the very strong mass dependence of these quantities, the heavy baryon 
result can be obtained for large nucleon masses, but the loop corrections change rather dramatically as the mass 
is varied. We mention in particular the complete resummation of the kinetic energy insertions in our approach, 
whereas these are only build up order by order in the heavy fermion scheme. Stated differently, the 1/m ex- 
pansion underlying HBCHPT is only slowly converging for the spin-dependent structure functions considered 
here. The delta Born and vector meson contributions to (Q 2 jmv) Si(y, Q 2 )|j, = o are displayed in Fig. [?]. First, 
these contributions lead to a visible change in particular for the proton and second, the uncertainty is smaller 
as for S\ '"(0, Q 2 ) chiefly because the second structure function does not depend on the parameter Y Q. Again, 
the inclusion of Air loops is called for. 

Recently, data from Jefferson Lab on the momentum dependence of the generalized DHG sum rule for the 
neutron at low and intermediate Q 2 have become available [[l2| . In Fig. || we show the lowest two data points 
for the integral /^(Q 2 ) in comparison to our theoretical prediction. Consider first the left panel, which refers to 
the results without the form factor of Eq.( |3.14 ). As stressed before, the theoretical uncertainty is rather large 



due to the large range of values of delta parameters. It is, however, remarkable that the trend of the data is 
well described for values of g2,X,Y at the upper edge of their allowed values (in magnitude). Positive values 
for X, Y and the smaller values of g-i are clearly excluded from this comparison. For illustration, we show that 
with the set of parameters 52 = 9.5, X = Y = —0.725 one obtains a good description of the data (see the 
dot-dashed line). In the right panel of that figure, we show the results if in addition the form factor is included. 
As expected, the allowed band shrinks and furthermore, for the second set of delta parameters the theoretical 
description turns out to be in rather good agreement with experiment for Q 2 < 0.15 GeV 2 . The bending of the 
data to lower values with increasing virtuality cannot be explained any more, which presumably points towards 
the inadequacy of the form factor Eq.(3.14) for analysing the new precision data.^] It is also interesting to 



compare the chiral loop prediction for the integrals Ia(Q 2 ) and Ic(Q 2 ), as shown in Fig. |9j for the neutron. 
While the third and fourth order contributions are rather different, the total one-loop prediction for I^(Q 2 ) 
and Iq(Q 2 ) comes out almost equal, as expected, see Rcf. |l6| . 

We now turn to the BC integrals hiQ 2 ) for the proton and the neutron as shown in Fig. [h]. First, we note that 
the validity of the BC sum rule was shown for HBCHPT in [|| , it is also fulfilled if one uses infrared regularization 
(which is simply a consequence of some general principles like e.g. crossing symmetry). In Fig. [l^ we show the 
result for I2 from the IR calculation of jl3) where the form factors are taken up to 0(q 4 ) (corresponding to the 
curve denoted as IR: full). The difference with the dispersion-theoretical result [|8| is mostly accounted for by 
vector mesons and one therefore obtains a rather good description of I2 in IR plus vector mesons. However, 
strictly speaking this does not correspond to an 0(q 4 ) calculation of I2, in which one only retains terms in the 
product of the form factors up-to-and-including fourth order (corresponding to the curve denoted as IR: 0(<7 4 )). 
It turns out that this 0(q 4 ) calculation is rather far from the empirical results since higher order terms are 



important. We come back to this issue in Section 4.2. In fact, it turns out that the HBCHPT result of Ref. |8| 
looks better though individually the form factors are described slightly worse. ^ However, in that paper the 



Dirac form factor is simply a constant (that is the charge of the nucleon) and the complete Q 2 -dependence is 
thus given by the Pauli form factor, quite in contrast to our calculation and the empirical result. Note also that 
the HBCHPT result shows no curvature in contrast to the empirical curve. Given these observations, we do not 
share the opinion of the authors of [|| that the BC sum rule is a good observable to bridge the gap between low 
and high Q 2 . 

In Fig. |ll], we show the forward spin polarizability for the proton and the neutron at finite photon virtuality. 
For the proton, the chiral loop contribution is positive and increasing with increasing Q 2 , while the delta Born 
graph contribution is negative and essentially flat. In total, the sum is positive, at odds with the measured value 

* 7 We remark that the theoretical band based on our preliminary calculation shown in Jlj| is somewhat broader because we 
allowed for the delta coupling 32 = X = Y = which we now consider excluded. 

# 8 Note that one has to be careful when comparing the large mass limit from the IR calculation with the HBCHPT results. For 
example, the LECs used in do not correspond to the ones used in the HBCHPT calculation of test] and thus a seemingly 
contradicting result emerges. 
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of the spin polarizability (at the photon point) |l4| . At this point we can only speculate that a combination of 
higher order and delta lo op e ffects will account for the discrepancy. This deserves further investigation (see also 
the discussion in Section 4.2). The chiral expansion of 7q will be discussed below. For the neutron, the results 
are markedly different. While the chiral loop contribution decreases with increasing photon virtuality, the delta 
Born contribution rises, leading to a slightly decreasing sum. Note also that 7q is negative in contrast to the 
positive 7q, our theoretical bounds are: 



1.0 • 1CT 4 fm 4 < 7 £ < 2.6 • 1(T 4 fm 4 
-1.8 ■ l(T 4 fm 4 < 7™ < -0.2 • 10~ 4 fm 4 



(4.1) 
(4.2) 



One should not be worried about the difference of these numbers to the results presented in j| because in that 
paper a third order dimensionally regularized relativistic calculation was performed, which has no consistent 
underlying power counting. It is interesting to consider the chiral expansion of pion loop contribution to the 
forward spin polarizability of the proton and the neutron. These can be given in closed analytical form, ^ 



p 
To 



p 2 n 2 
e 9 A 



1 - 



1 f 9 

- (21 + 3k v + k s ) tt^i + ( (20 + -(k v + K s )) n 2 



\ 15 

+ (26 + 7n v + 7n s ) In /ij - — (125 + 49k„ + 63k s ) tt/z 3 
4^45 - 8M + 6.03 + 3.22 +0(fj, 2 ) 

4.64 , 



(4.3) 



76 



96tt 3 A/2F2 



1 



- (9 + 3k v — k s ) 7r [i + ( (5 + 3k v — 3k s ) fi 2 
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+4 (l + k v — k s ) fi 2 hifij — + 15/tt, - 21k s ) 7t /x 3 
4A5 - ^25 + ^00 + 0^68 +0( fJ 2 ) 

1.82 , 



(4.4) 



where the last number refers to the full (not expanded) result and all numbers arc given in the canonical units of 
10~ 4 fm 4 . The first term in these series agrees, of course, with the result of j22|. More interesting arc, however, 
the higher order terms which have stirred some controversy within HBCHPT calculations, see J|(^, |3l], [32[ [34j 
(more precisely, this refers to the terms ~ fi _1 ). Let us consider first the proton. Upon inspection, the chiral 
expansion does not seem to converge. However, we note that the sum of the four first terms in the expansion is 
5.39 • 10 -4 fm 4 , the difference with the full result coming essentially from the 0(^ 2 ) term, which brings another 
—0.64- 10~ 4 fm 4 . Naively, one would expect that the expansion parameter is [i ~ 1/7, but clearly the appearance 
of large numerical prefactors, in particular due to the large value of k v ~ 3.7, invalidates such a simple picture. 
Things look somewhat different for the neutron. The sum of the first four terms amounts to 1.87 • 10 _4 fm 4 , 
only a few percent off the full result. Note, however, the almost complete cancellation between the first and the 
second term. Therefore, to make a reliable prediction for the chiral loop contributions to 7q , one has to account 
for the terms of order /i as done here. To get an idea about the theoretical uncertainty for these numbers, we 
have replaced the physical values of n v and k s in the above expressions by the values of the LECs eg and cj 
to third order as given in fll3| |. To the accuracy we are working, this is legitimate. The resulting values for 
the forward spin polarizabilities are 7q = 4.95 • 10 _4 fm 4 and 7° = 1.70 • 10~ 4 fm 4 , not very different from the 
results given above. We have also identified the reason for the bad convergence, it is entirely related to the 
contributions from the Born terms, as suggested by the method used in p2] |. We thus give here the third and 
fourth order pole and non-pole contributions to 7q (note that to third order one has no pole contribution in 

# 9 Note that the terms of order pP (and higher) will be modified by two-loop corrections (but not the terms ~ fi~ 2 and ~ M~ x ). 
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the case of the neutron): 



p. 3, pole 

To 

p, 3, non — pole 

To 

p. 4, pole 

To 

p,4,non — pole 

To 

n, 3, non — pole 

To 

n. 4, pole 
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n, 4, non — pole 
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-0.99 


+ 0.24 


-0.04 


-0.02+ . 


. =-1.03 



(4.5) 



(4.6) 



where the first, second, third, . . . term refers to the terms of order 1/M 2 , 1/M ff , M°, . . ., respectively, and the 
last one is again the full (unexpanded) result at the given order. All numbers arc given in canonical units. One 
clearly sees that the pole graphs are entirely responsible for the convergence problems, whereas the non-pole 
terms show a good convergence at third as well as at fourth order. One thus might entertain the possibility of 
rearranging the power counting accordingly. We come back to this point in Section |l.2j . Finally, note that a 
complete one-loop calculation to 0(q 4 ) within HBCHPT only gives the terms up-to-and-including C(/i _1 ). 

In Fig. [l^ the longitudinal-transverse spin polarizability for the proton and the neutron at finite virtuality 
is shown. Here, the delta Born graph contribution is small, as can be easily understood making use of the 
multipole decomposition of Ref. |Q. Therefore, the chiral prediction for ^'"(Q 2 ) is certainly more reliable 
than for the forward spin polarizabilities. For 8^{Q 2 ) as well as for 8q(Q 2 ) the slope at Q 2 — is negative, 
consistent with the HBCHPT result of B. At the photon point, we have: 



1.7 ■ 10~ 4 fm 4 < 81 < 2.2 • 10~ 4 fm 4 , 
2.3 ■ 10~ 4 fm 4 < SS < 2.8 • 10" 4 fm 4 . 



(4.7) 
(4.8) 



We note that these numbers are not very different from the leading 1/M 2 singularity, which gives 8q = 8q 
2.23 • 10 -4 fm 4 . This can be made more transparent by considering the chiral expansion of these quantities, 
find 
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^23 
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/! In /i 

- 0^ 
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15(k v + 4k s ) I it n' 



(4.10) 



where again the number in the last row refers to the full one-loop result and all numbers are given in units 
of 10" 4 fm 4 . The chiral expansion of 8q is very well behaved, furthermore, the sum of the first four terms in 
the series is 2.13 • 10~ 4 fm 4 , only 4% off the full result. For the neutron, we observe that the sum of the first 
four terms, which is 2.66 • 10~ 4 fm 4 , essentially amounts to the total sum at one loop. The chiral expansion 
of 8q is well behaved, although the terms of order /x° are somewhat small. We note that in contrast to the 
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chiral expansion of 7q : ™ no unnaturally large prefactors appear here. It thus appears that the longitudinal- 
transverse spin polarizability is a better suited observable to test the chiral dynamics of QCD. Finally, we have 
also calculated the longitudinal-transverse spin polarizabilities using the third order values for C6,7 instead of 
the physical values for k s ,v This results in Sq = 2.45 • 10~ 4 fm 4 and Sq = 3.13 • 10~ 4 fm 4 . 

4.2 Beyond fourth order 

As we have seen e.g. in the discussion of the forward spin polarizabilities, the calculation of the Born graphs 
to fourth order is not quite sufficient in some cases to obtain convergence. This observation was also behind 
the modified method to obtain the spin polarizabilities from the spin-structure amplitudes used in j32|. We 
therefore assume that the missing physics largely resides in these Born terms, more specifically in the leading 



two-loop corrections that appear at fifth order While a complete fifth order calculation goes beyond 
the scope of this paper, it is straightforward to calculate the contribution from the pole diagrams shown in 
Fig. |l3| (for the combination of structure functions S\{v, Q 2 ) — {Q 2 /mv) 52 (f, Q 2 ))- This class of diagrams 
is gauge-invariant and, of course, does not modify any LET or the DHG and BC sum rules. However, new 
LECs related to dimension two insertions as in diagram 5t and dimension four insertions as in diagram 5u 
appear. These have, however, already been determined in the analysis of pion-nucleon scattering |56| and of 
the nucleons' electromagnetic form factors jHS. We use the same values, i.e. C4 = 3.4 GeV -1 , e 54 = 0.26 GeV -3 
and 674 = 1.65 GeV~ 3 . Other operators that appear in the chiral expansion of the nucleon from factors [^3| 
do not contribute to the particular combination of the structure functions we are considering here. Note that 
here we concentrate on the neutron, simply because there are much less diagrams to calculate than in the case 
of the proton due to the non-vanishing charge (dimension one insertion). For further calculational simplicity, 
we restrict these considerations to observables that can be derived from the combination of structure functions 
Si(v, Q 2 ) — (Q 2 jmv) S%{v, Q 2 ). Of course, one should devise a modified power counting to properly account for 
these effects. A resummation of these terms also to higher order could be obtained by including the pertinent 
higher order insertions in the baryon propagator. Whether this can be done in a systematic fashion consistent 
with already existing successful calculations of other observables remains to be seen. Here, we are only interested 
to investigate whether these terms indeed contribute significantly in the problematic cases, these are the forward 
spin polarizability and the convergence of the BC sum rule. 

Consider first the BC sum rule for the neutron. Here, we have performed the complete fifth order calculation 
since only Born terms contribute. As shown in Fig. |lj, there are sizeable corrections when going from fourth to 
fifth order in the expansion of the product of the form factors. On the other hand, for Q 2 < 0.15 GeV 2 , the fifth 
order result is not very different from the "full fourth order" result (meaning that the form factors are taken to 
fourth order, not their product). Finally, we turn to the forward spin polarizability of the neutron. Using the 



same notation as in Eq.(4.6), we find (in canonical units): 

7o' 5 ' polc = +0 +0.02 -1.36 -0.18 -0.09 +... =-1.41, (4.11) 

where the first two zeros are consistent with the power counting, i.e. that the two leading chiral singularities 
~ 1/M 2 . and ~ are not affected at two (and higher) loop accuracy. Thus, the Born contribution gets 

sizeably reduced at this order. It is also interesting to note that the formally leading term ~ /i° is very small. 
Since the proton pole contributions have the same sign as the ones from the neutron and are larger in magnitude 
(related largely to the fact that the magnetic moment is larger), one can speculate that these terms will sizeably 



reduce the positive prediction, cf. Eq.(4.1). Such a calculation should be done but goes beyond the scope of this 
article. Finally, we have considered the moment I^iQ 2 )- Although we find some sizeable fifth order corrections 
from the Born terms, these are counterbalanced to a large extent by vector meson contributions, so that the 
ensuing band for the theoretical prediction is not very far from the fourth order one as given in Fig. 



# 10 For fifth order calculations within HBCHPT, see 35 
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A Chiral contributions to the structure functions 

Here, we give the explicit expressions for the reduced invariant functions A(s 1 Q 2 ) and B(s,Q 2 ) as defined in 
Eq.(2.7), for the third and fourth order contributions. Throughout, the contribution from the elastic intermedi- 
ate state is subtracted. For notational simplicity, we do not display the bar as done for the structure functions 
Sf'^ {v, Q 2 ) in the main text, cf. Eq. (2.6). Note that the B(s,Q 2 ) as given below correspond to the quantity 

"M", Q 2 ) = vS 2 {v,Q 2 ) - vS 2 (v, Q 2 )L =0) 

that is they allow to calculate the moments of S 2 for i > 1. However, from the 0(q 4 ) graphs which have an 
elastic piece, i.e 4(a)-4(j), one gets an additional inelastic contribution. Indeed, one has: 



B(s, Q 2 ) - B(2m 2 - 2Q 2 - s, Q 2 



-B(Q 2 



(A.l) 



where vb — Q 2 /2m and the ellipsis denotes the inelastic contribution given below. Using the fact that the pole 
contribution to vS 2 reads 



vS 2 {v,Q') 



2 \ pole 



2v 2 -v 2 B 



F 2 (Q 2 m(Q 2 )+F 2 (Q 2 )) 



(A.2) 



one gets: 



^- 1} (0,Q 2 ) - B(Q 2 ) = -F 2 (Q 2 )(^(Q 2 ) + F 2 (Q 2 )) 



(A.3) 



It is straightforward to check that this relation holds. 

Consider first the third order diagrams given in Fig. |l| and their crossed partners. We have (whenever appro- 
priate, we add together the contributions from various diagrams): 



3 a 



A{s,Q 2 ) 



2 2 
me ■ g% r 



2F 2 



1 

2m 2 



Ms) + 7T^(s + m 2 + Q 2 )Ji{s) 



B(s,Q 2 ) = 



(A.4) 



3b+3c 



A(s,Q 2 ) 



B(s,Q 2 ) 



2 2 
me g A r 



F 2 



Ms) 



1 



2 

2„2 r 



2m 2 



(s + m 2 + Q 2 )Ji(s) - 2Q 2 T 1 (s, Q 2 ) 



( s _ m 2 + Q 2 )T 2 (s, Q 2 ) + 4r 3 ( S , Q 2 )j -(1 - r 3 ) , 



2me g A 
~~F 2 



r 1 ( s ,Q 2 )-r 4 ( s ,Q 2 )l i(i-r 3 ) 



(A.5) 



3d+3e 



A(s,Q 2 ) 
B(s,Q 2 ) 



4me 2 g A 
F 2 " 



73(5, Q 2 ) 



2 2 

me z g A r 



F 2 



2 74 ( S ,Q 2 )+ 7 i(s,g 2 ) 



(A.6) 



3f 



A(s,Q 2 ) 



B(s,Q 2 



4:me 2 g\ 

F 2 

2 2 
me g A 

F 2 



-Ms, Q 2 ) + l(s-m 2 + Q 2 )n,(s, Q 2 )] i(l - r 3 ) 
2 74 ( S ,Q 2 )+ 7l ( S ,Q 2 ) + 2m 2 ^5(s,Q 2 )l - Vi) ■ 



(A.7) 
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3g 



A(s,< 



B(s,Q z ) = - 



me 2 g\ 




F 2 


■ 


+ (3m 2 - 


- Q 


-Q 2 (s- 


- to 


2me 2 g 2 A 


F 2 



| Jo 00 + i— a (* + m 2 + Q 2 ) Ji(s) + [s - m 2 + SQ^T^s, Q 2 ) 
! 4to z 

+ (3m 2 - Q 2 - s)T 2 (s, Q 2 ) - 4r 3 (s, Q 2 ) + (6 - d)(« - to 2 + Q 2 )G 3 (s, Q 2 ) 



r 1 ( s ,Q 2 )-r 4 ( s ,Q 2 ) 7 (3-r 3 ) 



(A.8) 



3h 
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3k+31 



3m+3n 



A{s,Q 2 ) 



A( Sl Q 2 ) 



me 2 g\ [ 3 



3Q 2 



F 2 
3Q 2 m 2 



A Ms) + — 2 ( S + m 2 + Q 2 )J 1 (s) 2(s _ m2) 
1 



(jo(s)-2J x ( S ) 



(s — m 2 y 



2 2 -1 

me g A r 1 



_F 2 



2TO 2 



(m 2 -Q 2 -s)Ji(s) 



S — 772 



A( S ,Q 2 ) = 



, , _ me 2 g\ ( 4Q 2 



F 2 



( J ^)7 3 (^Q 2 )J(l + -3) 
V s — m z ) — A 2 



(A.9) 



(Ji(s)-Jo(s))] ~(l + r 3 ) , (A.10) 



(A.ll) 



A( S ,Q 2 ) 



me 2 g\ r 1 



i?2 



, 2 (to 2 - Q 2 - s) Jx(s) - Q 2 r!( s , Q 2 ) + -(« - to 2 + Q 2 )r 2 (s, Q 2 ) + Q 2 r 5 (s, Q 2 ) 
4?n 2 2 



Q 2 



s — TO \ 2 



£ Ji(s) - ~ io(*) - (d ~ 2)E 3 ( S , Q 2 ) + Q 2 r 4 ( s , Q 2 ) 



2 

->2\ 2 



r 5 ( s ,Q 2 )-TO 2 r 6 ( s ,Q 2 ) -(i + r 3 ) 



i 



(A.12) 



3h+3i+3,j+3k+31+3m+3n 



B(s,Q 2 ) = 



me 2 g 2 A 
F 2 



7l ( S ,Q 2 ) + 2 74 ( S ,Q 2 ) + 



2m 2 



S — TO' 



-(7 2 (s,Q 2 ) + 27 5 ( S ,Q 2 ))] -(l + r 3 ) , (A.13) 



with 



F(s,Q 2 ) = F(s,Q 2 )-F(m 2 ,Q 2 ) , 
Z(s,Q 2 ) = Z(s,Q 2 )-(s-to 2 )F'(to 2 ,Q 2 ) , 



(A.14) 



where the prime denotes differentiation with respect to s with s = (p + q) 2 = 2mv + m 2 — Q 2 the standard 
Mandelstam variable and v is the energy transfer. The various loop functions are generalizations of the ones 
defined for Q 2 — in Ref. Jl^] to virtual photons with Q 2 > using 



h^x, y; s, Q 2 ) = M 2 (l - y) + M 2 y 2 + (s - M 2 )xy(y - 1) + Q 2 x(l - x)(l - yf 



(A.15) 



in their Feynman parameter representation. Note that the result for A(s, Q 2 ) agrees (up to an overall factor of 
m) with the relativistic one given in Appendix A of Ref. [Q if one uses the relations 

{s-m 2 + Q 2 )n 3 (s,Q 2 ) = 73 ( S ,Q 2 )+T 3 ( S ,Q 2 )-(s^2m 2 -2Q 2 -s) , 

2G 3 (s,Q 2 ) = T 4 ( Sl Q 2 )-T 2 ( s ,Q 2 ) , (A.16) 
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that hold in dimcnsionally regularized but not in infrared regularized baryon CHPT. 

Consider now the fourth order diagrams given in Fig. || and their crossed partners. After properly taking 
into account the renormalization of the anomalous magnetic moment fl3|| , we find (whenever appropriate, we 
add together the contributions from various diagrams. In particular, we have not shown the diagrams which 
arise from a trivial reordering of the two vertices, like e.g. for the contribution from graphs 4e)+4f). Such 
contributions are included in the results given.): 
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(A.17) 
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-2M 2 r (s, Q 2 ) + 4r 3 (s, Q 2 )] + 2Q 2 (j^s) - 4m 2 ri(«, Q 2 ) 
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Figures 




Figure 1: Fcynman diagrams contributing at third order. Solid, dashed and wiggly lines denote nucleons, pions 
and photons, in order. All insertions are from the leading order (dimension one) effective Lagrangian. Crossed 
diagrams are not shown. 
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Figure 2: Topologically incquivalcnt Fcynman diagrams contributing at fourth order. Solid, dashed and wiggly 
lines denote nucleons, pions and photons, in order. The black filled circle denotes an insertion from the dimension 
two effective Lagrangian. All other insertions are from the leading order (dimension one) effective Lagrangian. 
Crossed diagrams are not shown. 
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Figure 3: Predictions for the structure function Si(0, Q 2 ) with the elastic intermediate state subtracted in units 
of GeV~ 2 . The solid lines show the fourth order (one-loop) result whereas the bands given by the dashed lines 
refer to the one-loop plus resonance (delta Born graphs and vector meson) results. Left (right) panel: Proton 
(neutron). 
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Figure 4: Prediction for the moment T^Q 2 ). The solid lines show the fourth order (one-loop) result whereas 
the bands given by the dashed lines refer to the one loop plus resonance (delta Born graphs and vector meson) 
results. No NA'y form factor was used. 
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Figure 5: Chiral loop contribution to the first moment S^\o, Q 2 ) with the elastic contribution subtracted in 
units of GeV~ 4 . The solid (dashed/dot-dashed) line gives the result of the present calculation to order q 4 (the 
contribution from the third/forth order) in comparison to the heavy baryon result of || (dotted line). Left 
(right) panel: Proton (neutron). 
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Figure 6: Chiral loop contribution to the structure function (Q 2 /mi/)S2(0, Q 2 )|i/=o with the elastic contribution 
subtracted in units of GeV -2 . The solid (dashed/dot-dashed) line gives the result of the present calculation to 
order q 4 (the contribution from the third/forth order) in comparison to the heavy baryon result of (dotted 
line). Left (right) panel: Proton (neutron). 
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Figure 7: Chiral loop (solid line) and one-loop plus resonance contributions (dashed lines) to the structure 
function (Q 2 /miy) §2(0, Q 2 )\v=o with the elastic contribution subtracted in units of GeV -2 . Left (right) panel: 
Proton (neutron). 
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Figure 8: The integral Ia(Q 2 ) for the neutron in units of ^b. The solid line gives the fourth order result, 
the das hed lines represent the theoretical uncertainty due to variation in the delta parameters as explained in 
Section ^2. For comparison, the HBCHPT result is depicted by the dotted line. The data are from Ref. []l2|| . 
Note that here we have adjusted our normalization to the one used by the experimentalists. Left (right) panel: 
Without (with) AiV7 form factor. The dot-dashed line in the left panel corresponds to the delta parameters 
.92 = 9.5, X = Y = —0.725 as discussed in the text. 
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Figure 9: The integrals Ia(Q 2 ) — Ia{0) (left panel) and Ic{Q 2 ) — ^c(O) (right panel) for the neutron in units 
of GeV -2 . The solid line gives the fourth order result, the third and fourth order contributions are depicted by 
the dashed and dot-dashed lines, respectively. 
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Figure 10: The integral hiQ 2 ) for the proton (left panel) and the neutron (right panel). The solid line gives the 
complete fourth order result, that is the product of the form factors is taken to that order. Dashed (dot-dashed) 
line: taking into account the full fourth order (dispersion-theoretical) result for the form factors (i.e. without 
truncating the product of the form factors). The HBCHPT result of Q is also shown (dotted line). 
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Figure 11: The forward spin polarizability 70 at finite virtuality for the proton (left panel) and the neutron 
(right panel) in units of 10 -4 fm 4 . The solid line gives the fourth order pion loop result, the delta contribution 
is given by the band spanned by the dashed lines. The total prediction is depicted by the dot-dashed lines. 
The MAMI data point on j$(Q 2 = 0) is shown by the circle JTJ]. 




Figure 12: The longitudinal-transverse polarizability 80 at finite virtuality for the proton (left panel) and 
the neutron (right panel) in units of 1CP 4 fm 4 . The solid line gives the fourth order pion loop result, the 
delta contribution is given by the band spanned by the dashed lines. The total prediction is depicted by the 
dot-dashed lines. 
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Figure 13: Topologically different pole diagrams for V 2 CS off the neutron contributing at fifth order to the 
combination of structure functions Si(v, Q 2 ) — (Q 2 /mv) S2 (v, Q 2 ). Solid, dashed and wiggly lines denote nucle- 
ons, pions and photons, in order. The black filled circle (square) denotes an insertion from the dimension two 
(four) effective Lagrangian. All other insertions are from the leading order (dimension one) effective Lagrangian. 
Crossed diagrams are not shown. 
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Figure 14: The integral hiQ 2 ) for the neutron. The dashed lines represents the full fourth order results (that 
is the form factors taken to fourth order). The solid/dot-dashed/dotted line gives the fifth/fourth/ third order 
result (i.e. the product of the form factors taken to that order). 
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